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Direct Observations of Excess 
Solar Absorption by Clouds 

Peter Pilewskie* and Francisco P. J. Valero 

Aircraft measurements of solar flux in the cloudy tropical atmosphere reveal that solar 
absorption by clouds is anomalously large when compared to theoretical estimates. 
The ratio of cloud forcing at an altitude of 20 kilometers to that at the surface is 1 .58 
rather than 1.0, as predicted by models. These results were derived from a cloud 
radiation experiment in which Identical instrumentation was deployed on coordinated 
stacked aircraft. These findings indicate a significant difference between measure- 
ments and theory and imply that the interaction between clouds and solar radiation is 
poorly understood. 


iotis of abs<.)rption by pcrliaps ^5 W 
m o\'cr the Pacific warm pool (3). The 
relative error this difference introduces into 
current the(.>retical estimates of solar ;ib- 
sorptivu'i is lar^e, considering that avera^^e 
clear-sky absorptioi^ in tliai region is about 
100 W m ‘h The abstilute error appears to 
he sm<ill when compared to either terms in 
the etiermy bud^Tt, but that is misleading. 
Most of the solar radiation absorbed in the 
tropics ut^es toward heatini^ the surface; the 
remainder, about 20Vo, lielps dri\e the at- 
mospheric circulation. Thus, what a[^pear 
tt> be small errors in absorption by the at- 
mosphere mi^ht have hu^e consequences in 
tropical atmospheric dynamics. Another 
consequence of our inability to predict the 


evidence from several experimental and 
theoretical investigations <ner the past tour 
decades has shown that the ma<^nitude of 
shortwave (solar) absorpritm by clouds is 
uncertain. There has been some hitit that 
soLw absorjstion is in excess of that predict- 
ed by miKlels (1). C'ess et uL (2) and Ra- 
manatltan ct ai (3) reported tliat the ab- 
sorption by the entire atmospheric column 
in the presence of clouds exceeds model 
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mapnitude of solar absorption hy clouds is 
the misinterpretation of remote sensinL; 
data u.sed to infer c loud micTophysical prop- 
erties. in this report, we present measure- 
ments of cloud ahsorptiotr from the Tnqsi- 
cal CX'ean (jlohal Atmosphere ('oupled 
C\e<m Atmosphere Response Hxpenmeni 
(TOC jA-CX 3ARE) and the C'entral Equato- 
rial Pacific Experiment (CT:PEX), Eased (^n 
direct observations frcun aircraft. For con- 
sistency, we present our analysis in a man- 
ner similar to that used in (2, 3). 

In TOGA-CCMRE and OEPEX, 20 co- 
ordinated fli^dais were ma(.ie with i<.lenrical 
instnimenration above (at an altitude of 
20 km) and beneath (H to 12 km) cloud 
layers ro derermine cloud energetics (that 
is, flux divergence, ahsorptioi^, heatnu;, 
and so on). From TCOGA-C XTARh, G 
hours of useful solar ra<.liation data were 
acejuired durinp» well-coordinated flight 
segments (aircraft within C15°). (d:PhX 
provided an additional 18 hour> of 
well-coordinated solar flux measuremenrs. 

Hurin^ both TCX j A-( X^ARE atul C 'EPEX, 
the National Aeronautics and Space Ad- 
ministration (NASA) Pdl-2 aircraft (lew 
at {\early constant altitude near the tropo- 
pause, approximately 20 km; in TOC'jA- 
( X")ARE the NASA IXX8 flevs at mid- 
troposphere altitude, between 8 and 12 
km, and in ('EPEX, the mid-troposphere 
aircraft was the Aeromet Learjct. Each 
aircraft was instrumentc\l with two iden- 
tical broadband (CX^ to 4 0 |jun) sokir 
hemispheric field-of-view ra^licanclers 
(BBIIFC^VO fnr simultaneous measure- 
ment of upwcllin^ and downvudlinp flux 
at both flij^lir levels. Total direct-(.l ittuse 
ratliometers (TDDR) on each aircraft were 
used to measure spectral components of 
the sol ar flnx (4, 

It one is to determine the ahstirjstion in 
a layer, net solar flux must he acquircxl 
simultaneously, <ir nearly so, ar both flight 
altitudes (Fi^. 1). Usinj^ the flight navipa- 
tional data from the IX '-8 (or Learjet, for 
CdTEX) and ER-2 aircraft, wc shifted time 
series of flux data to best ali^n the data sets. 
Typical time offsets between ER-2 and 
IXT8 (Learjet) data were less than ^ min, 
ansi in most eases the offset was ne^liL^ihle. 
Therein lies the advance in the TCX^A- 
( XIARE and C'EPEX data sets over ths‘ data 
obtained earlier. Most of the earlier exper- 
imental attempts at determining cloud nh- 
>orption relies] on a single aircraft inakinp^ 
measurements ar several flight altitusles. 
The reslucrion ro elousl absorption then 
relicsl on knowleslj^e of cloud advestion, 
homogeneity, evoluri(»n, ansi so forth. Be- 
cause flux divergence is ohtainesl from the 
resislual of two relatively lart.;e numhers, the 
net fluxes, coeval me<tsurements are s rueial 
to limiting errors. Some limitesi atts’mpts 
have been made at flying stacked aircratl 
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(6), bur iJcnticiil instniments vvirb iclcnriciil 
calibration were nor used. 

We reduced our data in a manner similar 
ro rhar described in (2, d) to examine col- 
umn absorption in the presei'ice ot clouds. 
An important Llistincfion betweeri our mea- 
surements and tliose ot (2) is that our net 
tlux below cloud Livers is between <S to 12 
km abtne tlie surface rather than at the 
surface. We adopted two methods tor de- 
ducing: the effect of absorption by clouds (2, 
.■?): (.letermination of (i) the ratio ot cloud 
foreim: above the cloud to cloud forcin^^ 
beneath the cloud layer an«.l (ii) the slope of 
cloud reflectance versus transmission. 

(doud forcinp^ is defined as the difference 
between cloudy and clear-sky net flux. 
Cdoud forcinv beneath the cloud l.iyer, ( 
is equivalent to the difference Ivtween ab- 
soiption by the surface and the lower atmo- 
sphere under the clouds on the one hand 
and the absorption by the surface and the 
lower atmosphere under a clear column on 
the other, (doud torcin^: aK)ve the cloud 
layer (in our case, at the tropical tropo- 
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Fig.1. Net solar flux measured at 20 km from the 
NASA ER-2 (solid curve) and at 10 km from the 
NASA DC-8 (dashed curve) over the tropical west- 
ern Pacific during the TOGA-COARE flight on 31 
January 1993, 
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Fig. 2. Three-rriirtufe averages of cloud absorp- 
tion determined by the difference between net flux 
at 20 km and at 10 km, from TOGA-COARE and 
CEF^EX, The solid line indicates our estimate of 
clear-sky absorption. 


]Muse), CL], is the difference between ab- 
sorptioit by the cloudy total-column atmo- 
sphere and surface versus that by the clear 
rot<il-column atmosphere and surface. If the 
cloud forcing of the atmosphere column 
above the under-cloud level is denoted by 
then 

c:, - Ca + c; (1) 

and 

Cs/Ci = c:s/(Ca b (G) (2) 

The c[uantity is given by 

C'a ^ A(cloud) — A(cle.ir) (^) 

where A is the absorption by tlie column. 
The ratio ot to (L| is unity when 
A(cloud) is et|ual to A(elear); model pre- 
dictions place this rati<i around unity, and 
rarely does it exceed 1. 2 (3, 7), indicating a 
negligible difference in absorption be- 
tween the clear and cloudy cUmosphere. 
An important conclusion presented in (2, 
3) is that the ratio might be substantially 
higher, implying that a signific.int discrej'^- 
ancy exists between model-predictCLl 
cloud (c^r more precisely, cloudy column) 
absenption and cloud absorption derivcLl 
from measurements. 

The most critical element nee<.led to es- 
timate cloud forcing is the clear-sky net tlux 
because it cannot he determined concur- 
rently with a cloudy-sky net flux measure- 
ment. Rather than use satellite-derivcxl 
Ljuantities, we determined clear-sky net flux 
from direct measurements [as in (2)]. We 
constructed a scatter diagram of all net flux 
(F„J d3t;, versus the cosine of the solar 
zenith angle, By definition, the maximum 
values of for a particular value of solar 
angle occur under clear condititms. Howev- 
er, under broken cK>udy skies, downwelling 
solar flux can exceed the .same flux under 
clear conditions if .scattering from the sides 
of clouds becomes significant. With that in 
mind, we performed a linear regression anal- 
ysis of the upper envelope of points to de- 
termine the relation between (clear) and 
cos(^). The geographical region ctivered 
during TOGA-(X3ARE and C.EPFX was ex- 
tensive, between 140°E to I80°E and 15°S 


to 0°S. However, the entire region is over 
the ocean surface and is part of the quasi- 
homogeneous warm-pool region, where sea- 
surface temperatures are between 100 and 
101 K. The variability in humidity profiles 
above 10 km would not significantly affect 
clear-sky net flux over the region. We esti- 
mate the uncertainty of our clear-sky net 
flux to be ±5 W m •. 

(T)mparison of the clear-sky absoq^ri^m 
between the two flight levels, approximately 
10 to 20 km (solid line in Fig. 2), and 1-min 
averages of all cadumn (cloudy and clear) 
absorption obtained from the difference of 
F , measurements shows that, on average, 
cloutl absorption far exceeds clear-sky ab- 
.soqstion in the same layer. The average of 
the absorption was 165 W m maximum 
absorptit^n approaches 10‘/o of the solar con- 
stant. A few of the cloud absorption points 
actually fall under the clear-sky absorption 
values, possibly the result of side-scattering 
from clouds producing an anomalously high 
value of IX'-(S (Learjet) 

The dependence of absorption on :: is 
vlifficult to decipher from the 1-min a\er- 
aged ^.lata (Fig. 2). To examine this depen- 
dence, we computed the average clcnid ab- 
sorption in 0.1 -wide cos(^) bins. Thexse val- 
ues (Fig. 1) show <1 noticeable trend toward 
increasing absorption wath increasing co- 
sine, at a rate steeper than that for clear-sky 
absorptitm. 

The predicted model estimate of CT/(L| , 
uniiy, is not attained unless clear-sky and 
ck)udy-sky absi^rption are identical. (3ur 
data (Figs. 2 and 1) shenv substantial differ- 
ences between clear- and cloudy-sky ab- 
sorption. However, the model estimates are 
based on surface cloud forcing, whereas we 
tletermined cloud h)rcing beneath clouds 
but at 10 km. Caxsine-averaged values of 
L| A sltown in Fig. 4, were computcxl 
in a manner similar to the cosine-average^l 
absorption in Fig. 5. The average over all 
angles is 1.68. Using calculations of stdar 
transmission in the lowest 10-km layer (8), 
we estimate that this is equivalent to a ratio 
( T(.surface)/( ' , t)f 1 .5. Our findings are thus 
consistent with those reported in (2, 3). 

A more direct ctimparison with (2, 3) can 


Fig. 3 (left), Averages in 
0. 1 -wide cosine bins of the 
absorption data shown in 
Fig. 2, The solid line indi- 
cates our estimate of 
clear-sky absorption. Error 
bars denote the standard 
deviation of the sample, 
a measure of the range 
of variability in observed 
cloud absorption. Fig. 
4 (right). Same as Fig, 3, 
except that the ordinate is 
the ratio C> - /C^ . 
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Ix' made only by interring rlie surface flux ar 
the base of the column beneath the aircraft 
flight tmek. Usini>: an estimate of the averaj;]!e 
clear-sky column absorption from (d) and 
making the plausible assumption that our 
estimate of average cloudy column absorp- 
tion is somewhat independent of total col- 
umn depth (9), we can make such an infer- 
ence. Because of Eqs. 1 and 

( -s — C^l — A(cKxid) + A( clear) (4) 

C^ur estimates of A(cloud) and CJy are 165 
W m “ and -115 W m respectively. 
Adopting the estimate of A(clear) over the 
tropical warm pool from (3), 100 W m x 
we arrive at ,C>^,/;C^ j :■ — 1.58. We conclude 
that our estimates of cloud absorption are in 
considerable variance with model predic- 
tions (Fig, 5). 

A further consecjuence of using flight 
data versus surface data is the variability in 
planetary albedo along a flight track. All 
of the TOGA-CXOARE flights were over 
the sea surface, and st) variability in the 
reflectance at the altitude of the LXA8 was 
a direct measure of variability in cloud 
optical thickness and cloud cover beneath 
the aircraft. We therefore used our mea- 
surement of albedo from the 1X^8 to dis- 



Fig. 5. Summary of estimates of Cc^./ Cy, deter- 
mined by this study, {2, 3), and model calculations. 


criminate between clear and cloudy con- 
ditions beneath the cloud deck, adopting a 
threshold of 0.15. Althoiigh the albedo ot 
the sea surface is considerably less, around 
0.06 tor overhead sun, completely cle.ir 
conditions almost never occurred in that 
region of the tropics, particularly in tlie 
boundary layer, where scattered cumulus 
clouds w'ere prevalent. 

Se\ eral of the surface sites use^i by ( v‘ss cf 
till. (2) had available only flux radiometers 
that viewed the zenith; thus, measurements 
of surface net flux w'ere not possible. 1 low- 
ever, they were able to determine the rela- 
tion between satellite-derived albedo at the 
top t)t the atmospliere (TOA) and the ratio 
of instilation at the surface to insolation at 
the TOA— that is, reflectance versus trans- 
mission. (In our case, below'-cloud traivsmis- 
sion was measured ix)t at the surface but ai a 
flight altitude of 10 km.) The slope of the fit 
of their data is 0.6, whereas model pre- 
dictions are closer to —0.8. The lower mag- 
nitude slope for the measured data indicates 
that solar absorption in the column was 
greater than typically derived in modeb. 

For compari.son, w'e examined tlie rela- 
tion between cloud reflectance, measured .it 
the tropcqsause from the ER-2, and cloud 
transmission, measured beneath cloud lay- 
ers from the lXb8 (Learjet), flying between 
8 and 1 2 km, using the 3-min averaged 
integrated solar flux data. Filtering of se- 
lected data (Fig. 6) was done on the l^asis of 
albcxlo measured from the IX T8 and ,i 
threshedd i)f 0.15. Ft>r the 200 cases satisfy- 
ing that condition, the resulting slope of the 
regression line (Fig. 6) is —0.48. The vari- 
iibility of albedo measured from the 1X2-8 
albedo makes our analysis slightly different 
from that shown in (2), where stationary 
ground sites with fixed albedo were consid- 
ered. The tendency tow’ard higher surface 
albedo wanild increase the TOA albedo in 
the set of plants arising from minimal 
cloudy conditions (to the right in Fig. 6). 
The consequence of measuring the trans- 
mission from the IX b8 (as opposed to mea- 
suring transmission at the surface) has a 
similar effect. Transmission at 10 km ex- 


ceeds that measured at the surface, and this 
would also tend to decrease ihe slope. 

Although the slope of our regression fit 
is less than that found by CXxss ct ciL (2), it 
is consistent w'ith their general result^. Cail- 
culations of the atmospheric transmission 
between the surface and 10 km (8) show 
th.U the ground-based slope ot -0.6 woiilil 
fall to approximately —0.5 if surface trans- 
mission were replaced by transmission at 10 
km. On the basis of these results as well as 
our estimate of the ratio of cloud forcing 
beneath clouds to that at the tropopaiise, it 
is evident that our f inklings are completely 
consistent wath those of (2, .■?) and that 
measured cloud absorption exceeils theoret- 
ical estimates. 

Several factors that might contrilnite to 
the discrepancies between measured and 
calculated absorption by clouds liave been 
postulated (J ). Cdne is cloud inhomogeneity 
(models are typically run tt>r jslane-]^arallel 
homogeneous cloud layers). Filtering of our 
data allowcxl us to examine cloud iiK'rphol- 
ogy and inhomogeneity. We averagi-d our 
flux measurements over bmii^ peril kIs to 
siiuH)th higher frecjiiency features that 
might (Xherwise leail to difficulties in inter- 
[xeting data sets from the two aircraft plat- 
forms. The spatial resolution of tin. ER-2 
and IX b8 data are not identical, (H\ mg to 
differences in the ranges of clouds from 
both aircraft, thus necessitating the use ot 
time-averaged data. One measure of cloud 
homogeneity aloi'ig the flight track is the 
standard deviation of the upwxdling flux at 
ER-2 altitude (20 km) over the 3-mm pe- 
riod. We segregates] the slata shown in Fig. 
6 by filtering at 10 and 5'X> standarsi devi- 
ation in the upvvelling flux at 20 km (Fig. 
7). The slopes of the new regressioti litu's 
are —0,54 and —0.61 for the 10 atid 5% 
standard deviation threshokls, respecriwdy. 
Undoubtedly, the effect of filtering is to 
rcsluce the cluster of lower albedo ca>es and 
to increase the slope of our fit. Whether 
this can be construesl as indicating an in- 
crease in absorption wath inhotnogeneity 
remains uncertain. The statisricai signifi- 
catice of the filteresl data set is s|uestiotiable 
because only 25% ot all cases met the 5% 
standard deviation criterion. 

If cloutl morphology is the key to en- 
hanced cloud absorption, slo broken and 
variable clouds really increase column ab- 
.siirption (by changing [dioton path lengths) 
or is the effect only apparent and a measure 
of the inadequacy of current models to pre- 
dict tlie scattering and absorption of real 
clouds.^ eXher issues also remain unresolved. 
Ho clouds absorb anomalously, oi does 
there exist a different anomahuis abMirber, 
such as a carbon-based aerosol.^ Sjxctral 
data from the TCXiA-(X2ARE experiment 
suggest that this is not the case and that the 
significant absorption ami scatiermg fea- 


Fig. 6 (left). The ER-2 albe- 
do or refleotance, R, as a 
function of DC-8 transmis- 
sion. T. The solid line indi- 
cates the linear regression 
fit. The slope is -0.48. 
Fig. 7 (right). Same as Fig. 
6, except that two more re- 
gression curves are plotted, 
derived from data that were 
filtered to remove values 
arising when the standard 
deviation of the ER-2 net flux 
exceeded 10 and 5%. The 
slopes of the new curves are 
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0.54 and 0.61 , respectively. 
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tures in the tropical atmosphere are water- 
like. Further examination of this point must 
rely on more spectral measurements, partic- 
ularly in the near-infrared, where liquid 
water absorption and ice absorption become 
important. 
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